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Abstract 



■ We discuss the possible consistency of the recently discovered Type la supernovae at 

z > 1 with models in which dark energy is strongly coupled to a significant fraction of 
dark matter, and in which an (asymptotic) accelerated phase exists where dark matter 
and dark energy scale in the same way. Such a coupling has been suggested for a 
possible solution of the coincidence problem, and is also motivated by string cosmology 
models of "late time" dilaton interactions. Our analysis shows that, for coupled dark 
energy models, the recent data are still consistent with acceleration starting as early 
as at z = 3 (to within 90% c.l.), although at the price of a large "non- universality" 
of the dark energy coupling to difi^erent matter fields. Also, as opposed to uncoupled 
models which seem to prefer a "phantom" dark energy, we find that a large amount of 
coupled dark matter is compatible with present data only if the dark energy field has a 
conventional equation of state w > —1. 



There is increasing evidence that our Universe is presently in a state of cosmic acceler- 
ation [1] (or "late-time" inflation), i.e. that its energy is dominated by a component with 
negative (enough) pressure, dubbed "quintessence" , or dark energy. The energy density of 
such a cosmic component is, at present, comparable with that of the more conventional 
(pressureless) dark matter component, although their evolution equations can be widely 
different, in principle. 

In order to alleviate such a "coincidence" problem [2], without fine-tuning and/or ad 
hoc assumptions on the viscosity of the dark matter fluid [3], a direct (and strong enough) 
coupling between dark matter (or at least a significant part of it) and dark energy has 
been proposed [4]. Phenomenological models of this type have been widely studied [4, 
5, 6, 7], and shown to be compatible with various constraints following from supernovae 
observations, cosmological perturbation theory and structure formation. They are also 
theoretically motivated by the identification of the dark energy field with the string-theory 
dilaton [8, 9], and by the assumption that dilaton loop corrections saturate in the strong 
coupling regime [10]. 

A key feature of the coupled dark energy scenario is that it may allow a late -time attrac- 
tor describing an accelerated phase, characterized by a frozen ratio Pdark matter /Pdark energy 
This occurs, in models where a rolling scalar field plays the role of the cosmic dark energy, 
when the dark energy has an exponential potential and is exponentially coupled to scalar 
dark-matter fields (such an attractor property is lost, however, in the case of Yukawa-type 
interactions to fermionic dark matter [11]). For a successful scenario it is also required that 
a) the present coupling of the scalar to some exotic dark matter component is strong enough 
and approximately constant, and b) the present coupling to ordinary (baryonic) matter is 
weak enough, in order to avoid testable violations of the equivalence principle. 

In addition, the cosmic dark energy field has to be either massless or ultra-light, with a 
mass m ^ Hq, where Hq ~ 10~^^ eV is the present Hubble scale. For a scalar field gravita- 
tionally coupled to particles of mass M, on the other hand, there are quantum (radiative) 
corrections to the mass [12] of order m ~ gM(A/Mp), where q is the dimensionless cou- 
pling strength and A/Mp is the cut-off scale (typically, 1 TeV), in Planck units. This well 
known perturbative argument seems to suggest that both the coupling q to ordinary bary- 
onic matter (for which M ~ 1 Gev), and the value of M for the coupled, exotic dark matter 
components (for which the models assumes g ~ 1), need to be fine-tuned at extremely low 
values [13]. It should be stressed, however, that scalars with masses comparable to the (four 
dimensional) curvature of the Universe could be safe from radiative corrections, thanks to 
a recently proposed mechanism based on the AdS / CFT correspondence [14] . 

Quite independently from theoretical motivations, and from possible interpretations in 
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a string cosmology context, in this work we consider the basic scenario [4] in which dark 
energy is parametrized by a (canonically normahzed) scalar field (j), with exponential self- 
interaction energy of slope /x-\/273, 

= Foe"''^^. (1) 

Here /x > and Mp = (SttG)"^/^ is the reduced Planck mass. Neglecting radiation, as 
we are concerned with the late evolutionary phase of our Universe, we can describe the 
remaining gravitational sources as a cosmic fluid of non-relativistic, "dust" particles, dis- 
tinguishing however a "coupled" and "uncoupled" component, with energy densities pc and 
Pu, respectively. Uncoupled matter (for instance, baryons) satisfies the usual conservation 
equation 

Pu + ?>Hpu = 0, (2) 
while the coupled evolution of pc and (f) is described by the system of equations 

p^ + ZHpc-(5sl^^Pc = Q, (3) 

^ + 3i/^+|^+/3/|^p. = 0, (4) 

with /3>0. The total energy density is finally normalized according to the Einstein equation 

?>MlH^ = p^ + Pc + Pu, (5) 

where = 4''^ /2+V{<p). To make contact with previous work, we have used the notations of 
[4, 5], but we stress that the above equations also describe the asymptotic "freezing" phase 
of a "running dilaton" model [8], provided we identify the slope ^ a/2/3 and the coupling 
parameter /3-\/2/3 of the present model with the slope A and the dilaton charge q of [8] as 
follows: 

The above equations can be derived (for a conformally flat metric) from the (Einstein- 
frame) action 



/ 



d x^J—g 



where V'c and tpu are the coupled and uncoupled matter fields, respectively. The parameter 
f3 is thus related to the (homogeneous) scalar charge density, SSc/S(f), of coupled matter, as 

[2 Mp SS, 
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As shown in [4, 8], such a system of equations is characterized by an asymptotic accelerated 
regime in which the uncoupled component has redshifted away, p„ = 0, and all other terms 
of the Einstein equations have the same scaling behaviour: 

Fcx(^2o^/90ocpcocif2oca-3(i+'"), y; = (9) 

(a wider class of scalar field Lagrangians leading to similar attractor configurations has been 
studied in [9]). Within this regime p^/pc = const, and = p^/SMpH'^, Qc = Pc/3MpiJ^, 
are also constant. The evolution is accelerated, with constant acceleration parameter, 

d H q-1 2(3-11 , ^ 
= IH = = flO) 

provided 2(3>fi (or g > 1 in the notations of [8]). The luminosity-distance relation, for this 
regime, is obtained from the Hubble function 

H{z) = Ho{l + (11) 

(see Eq. (9)), where z{t) = aQ/a{t) — 1 is the red-shift parameter. The corresponding 
Hubble diagram is thus completely controlled by the dark energy parameter w, and its 
compatibility with previously known Type la supernovac (SNc la) has been studied in [5]. 

Motivated by the recent important discovery of many new SNc la at z > 1 [15], the aim 
of this paper is to re-discuss a possible fitting of present supernovac data in the context 
of the above model of coupled dark energy. Instead of assuming the asymptotic regime 
as an appropriate description of our present cosmological state, here we consider a "quasi- 
asymptotic" configuration, describing a cosmological state in the vicinity of the attractor, 
and for which the uncoupled component pu is small, but nonzero. The function H{z), 
and the associated luminosity-distance relation, will then be obtained by perturbing the 
"freezing" [8] (or "stationary" [4]) configuration (9), to the linear order. Such a perturbed 
configuration will contain the (critical) fraction of uncoupled matter, as a parameter 
measuring the typical (phase-space) "distance" from the asymptotic attractor with r2„ = 
(a different model with coupled and uncoupled fraction of dark matter components has been 
recently considered also in [7]). 

The uncoupled matter density scales in time faster than pc, namely pu ~ ~ (1 -|- 
zY, exactly like the standard dark matter component. Unlike in conventional models of 
uncoupled quintessence, however, the fraction of uncoupled matter is not to be identified 
with the total present fraction of non-relativistic (baryon-l-dark matter) components, and 
is not necessarily close to 0.3. In the case discussed here, the present value of $7^ may range 
from the minimum $7^ = ^baryons — 0-05 (we will refer to this as the "minimal model"), to 
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a maximum amount consistent with the clustered fraction of dust energy density ^Ira- As a 
reference value, we will assume Qm = 0.3. 

Such a different interpretation of Qu, as wc shall sec, has two important consequences. 
The first one, already stressed in [5], is that the accelerated regime may have a longer past 
extension, starting at early epochs characterized by z > 1, which are instead forbidden for 
the standard (i.e. uncoupled) quintessential models [15]. 

The second one, which probably represents the main result of this paper, is that a 
small enough value of is consistent with present data only for it; > —1. If we take, for 
instance, Q,u < 0.2, we obtain w > —1 at the 99% confidence level. More conventional 
analyses, assuming a fraction of uncoupled dark matter near to or larger than — 0.3, 
seem to favour instead w < —1 [16], although u; > — 1 is certainly still consistent with 
data. Should future data confirm these results, we would be left with the choice between 
"phantom" models of dark energy [17], uncoupled to dark matter but plagued by severe 
quantum instabilities [18] (see, however, [19]), and dilaton-like models of dark energy [4]- [9], 
non-minimally coupled to a significant fraction of dark matter. 

In order to perturb the dynamical system of coupled equations (2)-(5) it proves conve- 
nient to introduce the dimensionless variables 

2 - 2 _ ^ 2 _ Pu 

Obviusly, flc = 1 — — — u^. Denoting with a prime the derivative with respect to 
the evolution parameter N = ln(a/ao), the coupled system of equations, for the model 
presented in the previous section, can then be rewritten as [4, 20]: 



x' 



-^(l - x^ + y^) X - fiy^ + P{1 - a;2 - j/2 - u'^), 

^-i^'-y')- (13) 



3 

u' = -u I X 



.1 



Since there is complete symmetry with respect to simultaneous sign inversion of [i and ^, 
we will consider the case /x > only. It is also useful to write the equation for the derivative 
of R: one has, from Eq. (5), 

|^ = -^(3 + 3x2-3y2). (14) 

By equating the right hand sides of (13) to zero one obtains a system of algebraic 
equations whose solutions are characterized by constant fractional densities. Of particular 



4 



interest is the freezing (or stationary) configuration, given by 

3 V9 + 4/9/i + 4/32 

This point exists for ^ > and /x > (— /3 + ^/l8+~^)/2, is accelerated for n < 2(3, and is a 
global attractor [4]. On this solution, the properties (9)-(ll) are valid. 

The energy budget of our present Universe includes however an uncoupled component 
of non-relativistic matter at least as abundant as baryons, namely ^ 0, > 0.05, which 
measures our "distance" from the solution (15). In order to study the transient phase of 
approach to the attractor we will assume that is small enough to be compatible with a 
perturbative approximation, and we linearize (13) around (15), obtaining 

Sx' = i9xl-3yi-4Pxo-3)^-[3xoyo + 2yoif^ + P)]Sy-Pu^ 

Sy' = {iJ-yo + ^xoyo)dx + {3x1 - 9yl + 2ij.xo + ■ (16) 

In this approximation, we will assume that scales in time following the asymptotic 
behaviour (9), namely 

= nu {a/aof" = 0,e=^"'^. (17) 

We may thus regard as a first-order "external source", reducing (16) to an inhomogeneous 
system. The general solution is given by 

5x = ciie^i^ + ci2e^2^ + 

5y = C2ie'^^ + C22e''^ + A2U^ . (18) 

The coefficients Cij characterize the solution of the homogeneous system with v? = 0. 
Their values depend on the initial conditions, while the exponents £i and £i are found by 
direct substitution to be 

4{P + f,r 4(/3 + m)' 4(/3 + /x)2 +4(/3 + /x)- ^ ^ 

The square root of 

D = 324 - 32p^n - 63//2 _ Apfi{8n^ - 27) - 4^^(16^^ - 45) (20) 

is always imaginary in the region 2f3 > ii, for which acceleration occurs. Therefore, the 
"homogeneous part" of (18) describes at the linear level the oscillatory behavior of the 
system while converging to the attractor. Such oscillations clearly appear also at the level 
of exact numerical solutions of the coupled equations, see for instance [4, 8]. However, the 
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amplitudes Cij are always subdominat with respect to the inhomogeneous terms of Eq. (18), 
as shown by direct numerical integration of the linear system (16). 

In describing the approach to the attractor the most relevant role is thus played by the 
inhomogeneous contributions, proportional to u^. The coefficients Ai and A2 of Eq. (18) 
do not depend on the initial conditions, and are found to be 



8/3//3 + 2(9 + 8/32)//2 + 4/3^(2/^2 - 9) - 9(9 + 4/?2) ' 

, ^ /3[2M/9 + /^)-9]V9 + 4/3(/? + /x) 

^ 8/?M=^ + 2(9 + 8/32)^2 + 4^^(2/32 -9) -9(9 + 4/32)- ^ ^ 

By setting Cij = we can then compute, from Eq. (14), the Hubble parameter, and we 
obtain, to first order, 

2^ = -3(1 + w)- en^ixoAi - yo^2)e='"'^. (22) 

£1 

Direct integration leads then to the following approximate result for the perturbed Hubble 
parameter: 



rN 

^3wN 



Hi 



- enuixoAi - yoA2) / dNe^ 
Jo 

(1 - J^K)a-3(i+-) ^ f^^^-3l ^ (23) 



where we have introduce the "renormalized" fractional density of the uncoupled component, 
Qr, defined by 

,yoA2-xoAi 



nR = nu{'i- + e), e = 6- 

ow 

(24) 



3w 

9(9 + 2(3iJ.) 



3 • 



2/x2(9 + 8/32) + 4/3/i(2/32 - 9) - 9(9 + 4/32) + 8/3/x' 

Note that Qji/Qy^ is always larger than (but near to) unity, for realistic values of /3,/i [4]. In 
Fig. 1 we compare the result (23) for with a numerical integration of the exact equations 
(13) which includes the homogeneous part of the solution: the errors turn out to be very 
small, less than 3% in the range z < 6. 

The recent compilation of SNe la data can now be used to put constraints on the two 
parameters and w entering the expression for H(a) in Eq. (23). Such an expression is 
formally identical to that of a conventional, perfect fluid model with a prcssureless com- 
ponent, Qr, and a dark energy component, 1 — Or, with constant equation of state. We 
have performed a numerical analysis with the gold sample of data presented in [15], and 
we have reproduced the same results as in [15] for constant equation of state, but we now 
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Figure 1: Approximations to (assuming P = A, /j, = 3,Qb = 0.05, Qc = 0.3). The dotted 
curve is the numerical result, the full curve the linear fit. In the inset, the relative error. 

interpret as the total fraction of uncoupled matter rescaled by the (model-dependent) 
factor 1 + e. 

In coupled dark energy models, characterized by Eq. (9), the effective equation of state 
of the asymptotic regime always corresponds to a parameter w > —1. We thus begin the 
analysis by restricting ourselves to this region of parameter space. We also notice that, in 
all subsequent plots, the likelihood function has been integrated over a constant offset of 
the apparent magnitude that takes into account the uncertainty on the absolute calibration 
both of the SN magnitude and of the Hubble constant. All results are therefore independent 
of the present value of the Hubble parameter. 

In Fig. 2 we show the likelihood contours in the plane {w, 0,^). The value Qpi — = 
0.05 lies near the limit of, but inside, the contour corresponding to the 90% confidence level. 
This shows that an uncoupled percentage of dark matter density of the order of the baryon 
density is not ruled out by present SNe la data. Moreover, if some fraction of dark matter 
is also actually uncoupled, the fit is generally improved. However, Qu cannot be arbitrarily 
close to ilm, where Qm = + — 0.3 is the total non relativistic matter density that 
can be estimated through the mass in galaxy clusters. Moving toward the higher confidence 
regions of Fig. 2 requires indeed a lower and lower fraction of coupled dark matter, but the 
latter must be compensated with higher and higher values of the coupling P in order to get 
an equivalent amount of acceleration. 

This point is better illustrated in Fig. 3, where the confidence regions are plotted in the 
P-^u plane for = 0.27. For ci given value of the externcil parameter ^rn-! 

in fact, any 

point of the w — plane can be related to a point of the /3 — 0^ plane, using Eq. (24) and 




7 



0.4 




-1 -0.9 -0.8 -0.7 -0.6 -0.5 
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Figure 2: Contour plot obtained from the SNe la of [15], cutting to w > —1. Likelihood 
contours, from inside to ouside, are at the 68%, 90%, 95%, 99% confidence level. Qr is the 
effective fraction of uncoupled matter evolving as (see Eq. (23)). The line labelled by 
Qb marks the present baryon fraction. 



the relation 

(25) 

derived from = 1 ~ ~ ^O' valid near the attractor (we took the root that gives 
positive acceleration) . Values of Jim different from, but near to, 0.27, lead qualitatively to 
the same result. 

We can see from Fig. 3, where has been fixed at a reasonable value, that we 
cannot arbitrarily shrink to zero the value of without increasing f3. A coupling f3 much 
greater than unity, on the other hand, leads strong effective violations of the gravitational 
universality in the dark matter "subsector" corresponding to ^Ic- Such violations are in 
principle constrained by present observations comparing the distribution of dark matter 
and baryons in galaxies and clusters [21]. The corresponding upper bounds derived on /3, 
however, depend on several assumptions concerning the dark matter distribution and are 
of limited generality. Moreover, the bounds in [21] are derived assuming that all the dark 
matter is coupled, while here we are taking the more general point of view that Qc may 
differ from Q^n- 

For all these reasons, a firm quantitative limit on /? is hard to be derived from present 
data, and will be discussed in a forthcoming paper [22]. In this paper we will not exclude 
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2 4 6 8 10 12 14 

Figure 3: The SNe la results projected on the plane (3,i^u, at fixed = 0.27. Likelihood 
contours, from inside to outside, are at the 68%, 90%, 95%, 99% confidence level. The 
dotted curves mark the parameter values for which Zacc equals 1, 2 and 3, respectively. 
Values of ^acc above 1 are within one a, and even z^cc = 3 is not excluded at more than 
90% C.I.. 

a priori large values of /?, taking into account, however, that the large couplings required 
for a good fitting to SN data could impose severe constraints on coupled models of dark 
energy. 

In Fig. 3 we have also plotted the curves Zacc marking the beginning of the accelerated 
regime, which are obtained from Eq. (23) as 

z^ci^R, w) = [{3w + l){nR - _ (26) 

It follows that even Zacc = 3 cannot be excluded at more than 90% confidence level, while a 
value of Zacc > 1 is well within the 68% confidence level. This shows that the claim of [15] 
(see, however, [23]) about experimental detection of deceleration at z > 0.5 does not apply 
to models of coupled dark energy. 

On a more qualitative level, we show in Fig. 4 the best fit of SNe la for three classes of 
models, compared to the dataset of [15]: the standard ACDM scenario, a varying equation 
of state with w = wq + wiz, and our coupled model. The plotted variable is the distance 
modulus A(m — M) = (m — M)modci — {m — M)Miine) as customary (see for instance [24]), 
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Figure 4: Best fits of the distance modulus A(m — M) for three different models: ACDM 
{Qm = 0.27, r^A = 0.73), coupled dark energy {w = —0.6, Or = 0.05 ) and a varying 
equation of state w = wq + wiz, with wq = —1.31, wi = —1.48, fixing Q^n = 0.27. The data 
represent the gold sample of [15]. 

where apparent and absolute magnitude are related by 

m-M = 51ogiodL(^), dLiz) = {l + z) r dz'H~\z'). (27) 

Jo 

As it can be seen, the coupled model is almost degenerate (difference less than Am = 0.1) 
with ACDM up to z ~ 1.5, and with the varying- u; model up to z ~ 1. However, in 
these three models, the acceleration begins at very different epochs {zacc = 0.75 for ACDM, 
Zacc = 0.69 for w = Wo + wiz, to be compared to Zacc = 3.53 for the coupled model). 

We can also perform the fit by including the region of parameter space with w < —1. 
In that case, the best fit moves to larger values of Qu (see Fig. 5). As a consequence, the 
baryon value Or ~ Oj, ~ 0.05 moves near the 97% confidence level. This shows, on one 
hand, that additional SNe la data may have the potential to disprove models in which all 
(or most of) dark matter is coupled to dark energy. On the other hand, this also shows that 
if the uncoupled fraction of dark matter is small enough, namely 0„ < Or < 0.2, then we 
are lead to the region of parameter space with w > —1, excluding "phantom" or "k-essence" 
models of dark energy [17], possibly associated to embarrassing cosmic violations of the null 
energy condition. 

Finally, we can also derive constraints on w alone by marginalizing over Or (i.e. in- 
tegrating the full likelihood over the domain of Or, see Fig. 6). In the minimal case in 
which the totality of (renormalized) uncoupled matter coincides with the baryon fraction, 
Ofe, we marginalize over Or = 0.05 it 0.01 (we disregard here the possible small difference 
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Figure 5: Contour plot obtained from the SNe la of [15], including w < —1. Likelihood 
contours, from inside to ouside, are at the 95% and 99% confidence level. The variables 
w, ilji, ilfe are the same as in Fig. 2. 

between r2„ and The best fit turns out to be = —0.6 ± 0.05, rather independently 

of the exact value of Qr. If we keep instead the range = 0.04 — 0.3, i.e. the full possible 
range of the uncoupled matter fraction, we then obtain w = —0.8 it 0.1. This shows that, 
with some finite fraction of coupled matter, the allowed range of values of w remains firmly 
anchored to the w > —1 region of parameter space, in contrast with the standard case in 
which Qr ~ 0.3ib0.05, for which the peak of the likelihood is at tt; ~ —1.1 and the likelihood 
content of the w < —1 region is 80%. 

In conclusion, we have shown that from the analysis of recent SNe la there is still some 
room for strongly coupled dark energy scenarios which can ease the coincidence problem, 
and are compatible with an earlier beginning of the accelerated regime. Although more 
conventional models fit better the data, they seem to point toward a w < —1 equation 
of state for the cosmic dark energy component. If this tendency is confirmed by future 
observations we will face the problem of a theoretical interpretation of such a strong negative 
pressure. "Phantom" models, in fact, seem to be plagued with severe quantum instabilities. 
Coupled dark energy models, in this respect, could provide an interesting alternative. 

We found that the "minimal model" in which only the baryons are uncoupled fits the 
SNe la at 2a if we restrict to w > —1, and at roughly 98% c.l. if we include the region 
w < —1. In both cases the best fit for the equation of state is w = —0.6 it 0.1. The 
universe may start the present acceleration at a very early epoch, z ~ 3, much earlier than 
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Figure 6: Likelihood for w marginalized on Q^, with priors in three different intervals. Full 
curve: 17/?^ G (0.04,0.3), flat prior. Dotted curve: 0,r £ (0.04,0.06), flat prior. Dashed 
curve: £ 0.3 it 0.05, Gaussian prior. The peak moves to w > —1 when is smaller 
than 0.3. The horizontal lines mark the 68% confidence interval. 

any uncoupled model. If some fraction of dark matter is let to be uncoupled as well, then 
the fit improves steadily as the fraction increases, but at the price of a large coupling /3, 
which in turn induces potentially harmful violations of universality at the level of large 
scale gravitational interactions. A modest increase in the data statistics could definitely 
reject the minimal model of coupled dark energy, thus forcing an admixture of coupled and 
uncoupled dark matter as an acceptable cosmological scenario. 
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